Excitations of strongly correlated lattice polaritons 
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We present an analytic slave-boson approach to calculate the elementary excitations of the Jaynes- 
Cummings-Hubbard Model (JCHM) describing strongly correlated polaritons on a lattice in various 
quantum optical systems. In the superfluid phase near the Mott transition we find a gapless, linear 
Goldstone mode and a gapped Amplitude mode corresponding to phase and density fluctuations, 
respectively. The sound velocity of the Goldstone mode develops a peculiar anomaly as a function 
of detuning at low densities, which persists into the weakly interacting regime of a polariton BEG 



The realization of a BEC of weakly interacting polari- 
tons in high-Q microcavities [![, i.e., quasiparticles which 
form when photons strongly interact with quantum well 
excitons, has triggered an immense interest in quantum 
condensed and coherent light-matter systems. The ad- 
dressibility of a cavity with external lasers, device inte- 
gration, and high tunability make these optical systems 
ideally suited for quantum information and simulation 
tasks. Furthermore, the inherently composite nature of 
polaritons promises rich and interesting new physics. 

One of the most exciting questions in this emerging 
field is whether one can realize a Mott insulator - su- 
perfluid (MI-SF) transition of strongly correlated polari- 
tons. The Jaynes-Cummings-Hubbard Model (JCHM) 
has been introduced to describe such a quantum phase 
transition of light in an array of coupled QED cavities, 
each containing a single photonic mode interacting with 
a two- level system (qubit) [2j. Meanwhile the JCHM has 
become a paradigm lattice model for polaritons describ- 
ing various systems at the interface of condensed matter 
and quantum optics Q . Proposals based on cavity QED 
arrays can be realized using atoms jj|, excitons [5j, or 
superconducting qubits [f|. In trapped ion systems, e.g., 
in a linear Paul trap, a transverse phonon assumes the 
role of the photon and couples to hyperfine levels of the 
ion via an external laser Q ; tunneling of phonons is me- 
diated by Coulomb interactions between the ions. Ultra- 
cold atom systems with two optical lattices that couple to 
the photon mode of a single QED cavity may also realize 
a JCHM Q. 

Theoretically, the phase diagram and the excitations 
in the Mott phase have been calculated accurately in 
D = 1,2, and 3 dimensions However, very lit- 

tle is known about the nature of the superfluid phase, 
except the existence of a Goldstone mode predicted by 
heavy Monte-Carlo simulations in D = 1 [14J . In this let- 
ter we present an analytic slave-boson approach valid in 
both phases. In the Mott phase we find gapped particle 
and hole modes in quantitative agreement with previous 
results obtained from a diagrammatic linked-cluster ex- 
pansion 12|. In the superfluid phase we find a gapless 



In addition, our analytic approach is ideally suited to 
discuss the behavior of excitations in the superfluid phase 
as a function of detuning, which provides an important 
experimentally accessible parameter. Detuning can be 
used to change (i) the nature of polaritonic excitations, 
i.e., the relative weight of its bosonic and atomic parts, 
and (ii) the strength of the effective repulsion between 
polaritons. We find that the phase diagram as well as the 
sound velocity of the Goldstone mode show a surprising 
anomaly when this parameter is varied: For high polari- 
ton densities, the size of the Mott lobes and the sound 
velocity decrease for any finite detuning, while at low 
densities they steadily increase when tuning through the 
resonance. Comparison to an effective Bogoliubov the- 
ory for the JCHM shows that this anomaly prevails into 
the regime of a weakly interacting polariton BEC [TtJ • In 
this regime, superfluidity has been demonstrated recently 
by measuring the sound velocity of exciton-polaritons in 
a single micro-cavity using angle-resolved photon spec- 
troscopy fljj ]. Similar experimental techniques should 
also be applicable to lattice models as discussed here. 

The Hamiltonian of the JCHM is given by 
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linear Goldstone mode and a gapped Amplitude mode 
corresponding to phase and density fluctuations (similar 
excitations were found in ultracold gases [Hj]). 



where hj c denotes the local Jaynes-Cummings Hamilto- 
nian hj c = u c a\ai+w x af + g(af ai + &i a\) with site 
index i, boson creation (annihilation) operators a\ and 
qubit raising (lowering) operators ' . The bosonic 
mode frequency is w c , the two qubit levels are separated 
by the energy lo x and the coupling is given by g (we set 
h = 1). We also assume that the total number of excita- 
tions, i.e., polaritons N = X^i( a I a « + cr i +<T i r ): * s conserved 
and fixed by the chemical potential /i. We note that 
even under more general conditions including the effects 
of an external drive and dissipation the underlying equi- 
librium model still captures the essential physics of the 
MI-SF transition as recently shown for the Bose-Hubbard 
model (BHM) 0. The third term in {!]) describes the 
derealization of bosons over the whole lattice due to hop- 
ping between nearest neighbour sites with amplitude J. 
It competes with an effective on-site repulsion between 
bosons mediated by the coupling g. This competition 



leads to Mott lobes in the quantum phase diagram Q. 

The on-site eigenstates of the Jaynes-Cummings 
Hamiltonian hj are labelled by the polariton number 
n and upper /lower branch index a = ±. The mixed bo- 
son (n, n — 1) - qubit (g, e) states define upper and lower 
polariton states 

\n+) = sin0„|n,#) +cos0„|(n- l),e) , 

\n-) = cos6 n \n,g) - sm6 n \(n- 1) ,e) , (2) 

with the mixing angle tan#„ = 2gy/n/(S + 2\ n ), Xn = 
y g 2 n + S 2 /4 and the detuning parameter <5 = uj c — lo x . 
The corresponding eigenvalues are 

< = ~(n - Uc)n - 6/2 + axn , 



(3) 
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|0 ,g) is a special 



The zero polariton state |0) 
case with eo = = 0. 

A convenient starting point for our slave-boson ap- 
proach is the polariton representation [l3| of the boson 
operator at = Y^ na v /n"^n-i^i m terms °f standard 
algebra operators P[J = \na)u(Q\ and matrix elements 
/£" = { n -l v\a\n a) with f° v = (y/^ + av ^/n^T) /2 
for ?i > 1 (f£~ = l/y/2) at zero detuning (6 = 0). In 
this new basis the JCHM becomes 
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The polariton operators obey bosonic commutation rela- 
tions if the constraint 
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is f ullfillcd at each site i. In [12| we have shown that the 
presence of the upper polariton branch (cr = +) leads 
to additional high energy conversion modes in the Mott 
phase with small spectral weight and bandwidth; we ne- 
glect the upper branch as well as particle conversion tun- 
neling from now on and drop the branch index cr. In order 
to calculate the phase boundary and static observables in 
the superfluid phase near a Mott lobe with filling n, we 
restrict the Hilbert space to states with n and nil bosons 
and make a Gutzwiller Ansatz for the ground-state wave 
function 

\i/>) =l[[cos(e)Pj +sin(9)( S in(x)Pli + cos(x)Pl)] |0) 
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where we also dropped the index n, i.e., P ia . . ;„_,_,,■ 
( a n+a = e a , and f^ a = f a . This yields the variational 
energy e va r = (ip\H\ip), 



e cos((9) 2 + sin(6») 2 (e_i sin(x) 2 + ei cos(x) 2 ) 



Jz/4 sin(2#) 2 (/o cos(x) + /-x sin(x)X 




FIG. 1: Phase diagram for the JCHM displaying the low- 
est three Mott lobes with n = 1,2,3. Dotted lines represent 
the critical hopping strength's J c /g, where chemical potential 
and detuning are chosen such as to fullfill particle-hole sym- 
metry. A finite detuning \S\ > decreases the critical hopping 
strength J c /g for n > 1, but the lowest Mott lobe (n = 1) 
steadily increases when tuning through the resonance (<5 = 0). 

which has to be minimized with respect to the variational 
parameters 9 and x- The lobe boundaries as shown in 
Fig. [1] are determined by the vanishing of the order pa- 
rameter 4> c = (ip\a\ip) = sm(6>) [f sin(x) + /i cos(x)] 2 /2. 
Our results for the quantum phase diagram agree exactly 
with those in [III [ill]- Here, we point out that the size of 
all Mott lobes with filling factor n > 1 decrease for any 
finite detuning \6\ > 0, while the size of the lowest Mott 
lobe (n = 1) increases steadily as the system is tuned 
through the resonance ((5 = 0). We explain this special 
behavior when discussing the nature of excitations close 
to the lobes below. 

In order to find the excitations we follow the proce- 
dure outlined in [l6[ for the Bose-Hubbard model and 
define a new set of operators tV = (G\,E\ i ,E\ j ) T , 
which is obtained from the original polariton basis = 
(Pj , -Pj_i, Pii) T via a unitary transformation Rt = TP^ 
with 



cos(0) sin(0) cos(x) sin(0) sin(x) 
T = | sin(0) cos(#) cos(x) cos(6>) sin(x) 
- sin(x) cos(x) 



(7) 



The operator G T creates a new vacuum state, i.e., the 
mean- field ground state \ip) = ILC^O), and E^E^ 
are orthogonal operators creating excitations above the 
ground-state. We express the Hamiltonian in terms 
of these new operators and eliminate d by using the 
constraint §5§ in the restricted Hilbert space Gi ~ 

1 — E^Eu — E\ i E2i- Expanding the square root ev- 
erywhere in the Hamiltonian to quadratic order in 
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yields, after a Fourier transformation, an effec- 

t 



tive quadratic Hamiltonian H e g = e va r + X)k ^k ^cff.k E k 



(6) where E = (£ lk , £ 2k , Sj_ k , ^-k) T and h e s,k is a 4 x 4 



3 




FIG. 2: Spectral gaps and sound velocity as a function of 
tunneling strength J/J c o at zero detuning 8/g = and for 
(a) p, = jU c o where /i c o denotes the critical chemical po- 
tential at the tip of the lobe with critical hopping strength 
Jc = JcO (top figure) and (b) away from the tip at p, = 1.2p c o 
with J c = 0.566 Jco (bottom figure). Shown are the gaps of 
particle (dashed) and hole (solid) modes in the Mott phase 
(J < J c ) and the gaps of the Amplitude mode (dashed) and 
the sound velocity of the Goldstone mode (solid) in the super- 
fluid phase (J > J c ). The insets show the corresponding ex- 
citation spectra at i) J = 0.5J c o ii) J = Jco hi) J = 1.5J c q iv) 
J = 0.566 Jco- At the phase boundary, the particle and hole 
mode of the Mott phase are identical with the Goldstone and 
Amplitude modes of the superfluid phase. At the tip of the 
lobe ii), where the polariton density can remain constant dur- 
ing the superfiuid-insulator transition, the Amplitude mode 
becomes gapless and linear (its mass vanishes). The sound 
velocity of the Goldstone mode remains non-zero, confirming 
a special point in the phase diagram with dynamical critical 
exponent 2=1. Away from the tip iv), the Amplitude mode 
remains gapped and the Goldstone mode becomes quadratic 
with a vanishing sound velocity corresponding to a generic 
dynamical critical exponent z — 2. This confirms the change 
of the universality class along the phase boundary found in 
E2;0! but put into question in 

matrix. The sum over k runs over the first Brioullin 
zone. The effective Hamiltonian can be diagonalized by 
a bosonic Bogoliubov transformation [2(| yielding 

^var T ^fluct \ 

]>>a(k)4 kC U (8) 

a=l,2 k 

with efl uc t a fluctuation-generated correction of the 
ground-state energy and creating excitations with 
energy e Q (k). The analytic expressions for the two modes 
e a (k) are rather lengthy and will be presented elsewhere. 
They depend on the variational parameters 8 and \- I n 
order to calculate the excitations at different points in 



the phase diagram, we first determine the optimal pa- 
rameters 9 op t and Xopt from a minimization of (|6]) and 
insert them into the corresponding algebraic expressions 
for the spectra e Q (k). The results are summarized in 
Fig. rj] In the Mott phase, where 8 opt — 0, we find two 
gapped modes corresponding to particle and hole excita- 
tions of lower polaritons. Our results agree exactly with 
the ones obtained from a strong-coupling RPA [12|. In 
the superfluid phase we obtain a gapless Goldstone mode 
which is linear for small k describing the propagation of 
phonons with a characteristic velocity c s and a gapped 
Amplitude mode, which describes local density fluctua- 
tions between normal fluid and condensate. Its existence 
has been predicted for cold atoms in an optical lattice as 
described by the BHM 0. 

In cavity QED systems temporal and spatial correla- 
tion functions are experimentally accessible and provide 
an important tool to characterize the nature of different 
phases. In particular, the sound mode is a key signature 
of superfluidity and an important quantity which signals 
the presence of interactions in a condensate. In Fig. |31 
we study the dependence of the sound velocity c s and the 
condensate density p c = \4> c \ 2 on detuning 8/g. 

In the strong-coupling regime (J ~ J c ) the sound ve- 
locity develops an anomaly vs. detuning at low polari- 
ton densities p = (N) /N s (with N s the number of lat- 
tice sites). While for p — 2 (and higher) c s decreases 
for any finite detuning \S\ > 0, it increases steadily 
for p — 1 when tuning through the resonance. To ex- 
plain this behavior, we define the particle-hole gap in the 
atomic limit U n (6) = e„+i — 2e n + e„_i. This provides 
a useful measure for the effective repulsive interaction 
between polaritons and its dependence on detuning near 
the Mott lobes, where the mixing of states with differ- 
ent polariton numbers is small. For the Bose-Hubbard 
model U n (S) = U is simply the Hubbard interaction 
parameter but for the JCHM the particle-hole gap de- 
pends on filling n and detuning 8/g. For n 3> 1 we find 
U n (S) - U n (0) = -(3/8n^)(S 2 /g) + 0{8 4 /g\ 1/n 7 / 2 ) 
and thus U n decreases for any finite detuning \S\ > 0, 
consistent with a decrease of the sound velocity. How- 
ever, for n = 1 we find U^S) - C/i(0) = 5/2 + 0{8 2 /g), 
i.e., the effective repulsion increases steadily when tuning 
through the resonance. This causes the anomaly in c s at 
low densities. The different behavior of U\ is rooted in 
the special nature of the zero-polariton state with eq = 
independent of detuning. This special behavior is also 
the reason for the decrease of the condensate fraction 
p c /p for p = 1 (see inset in Fig. [3]) and the increase 
of the size of the lowest Mott lobe (see phase diagram 
in Fig. [J). At higher densities the condensate fraction 
increases while the size of the lobes decrease for any fi- 
nite detuning, consistent with a decrease of the effective 
repulsion. We conclude that in the strongly-correlated 
regime both sound velocity c s and condensate fraction 
Pc/ P are dominated by the repulsive interaction between 
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12-10 
6/g 8/g 
FIG. 3: Sound velocity as a function of detuning S/g for fixed 
polariton density p — 1 (left figure) and p = 2 (right figure) 
near the Mott lobes with J = J c o(S) + 0.2J c o(0) (solid lines) 
obtained from slave-boson theory and deep inside the super- 
fluid phase with J = J c o(S) + 5J c o(0) (dashed lines) obtained 
from an effective Bogoliubov theory. Here, J c o(S) denotes the 
critical hopping strength at the tip of the lobe. The insets 
show the condensate fraction p c /p with the same notation as 
above. 

polaritons when S/g is varied. 

In the weakly interacting regime (J 3> J c ), the slave- 
boson approach is no longer a good approximation since 
states with different polariton numbers mix strongly. We 
use an effective Bogoliubov theory instead, originally de- 
rived for a generalized Dicke model describing localized 
excitons coupled to photons conhned in a single mi- 
crocavity with a gapped, quadratic photon dispersion 
e(k) = ujq + k 2 /2m with ujq = irc/L and m = ir/(cL) 
(L being the width of the cavity) [l7|. The JCHM maps 
onto the generalized Dicke model after a Fourier transfor- 
mation of the bosonic field operator and expanding the 
lattice dispersion for small k- vectors as e(fc) « Jz + Jk 2 
(for a hypercubic lattice with lattice constant a = I). 
Using the results in (l7j with the formal replacements 
loq = lj c — zJ and m = 1/(2 J), we can compare our find- 
ings in the strongly interacting regime with those valid 
for weak interactions shown as dashed lines in Fig. [3] We 
observe that the sound anomaly is still present deep in- 
side the supcrfluid phase, but slightly suppressed. This 
is due to the admixture of states with higher polariton 
numbers n > 1 for p = 1. On the other hand, the max- 
imum of the sound velocity at higher densities p = 2 is 
slightly shifted to positive detuning due to a weak admix- 
ture of polariton states with n — 0,1. On the contrary, 
the condensate density behaves different in the weakly 
interacting regime far away from the Mott lobes. It de- 
creases independent of the polariton density when tuning 
through the resonance and is thus no longer dominated 
by the effective repulsion but rather by the internal struc- 
ture of the polaritons. To understand this we have to con- 
sider the boson density pb = p n + Pc, which consists of a 
condensate p c and a normal part p n (depletion). In the 
weakly interacting regime, the depletion p n is small and 
the condensate density p c behaves similar to the density 



of bosonic excitations pb; both decrease when the en- 
ergy cost of creating a boson is increased with respect to 
the energy cost of an atomic excitation, i.e., when tuning 
through the resonance. 

In summary, the slave-boson theory for strongly corre- 
lated polaritons presented in this paper predicts the exis- 
tence of gapped particle/hole modes in the Mott regime 
and, besides a gapless, linear Goldstone mode, the ex- 
istence of a gapped Amplitude mode in the supcrfluid 
phase. The anomaly of the sound velocity of the Gold- 
stone mode as a function of detuning and polariton den- 
sity is rooted in the special nature of the zero polariton 
state within the Jaynes-Cummings ladder. It is present 
in the strong as well as weak coupling regime and may 
constitute an interesting experimental signature for the 
composite nature of supcrfluid lattice polaritons in future 
experiments with coupled cavities. 
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